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Abstract. We report bulk superconductivity induced by an isovalent doping of 
phosphorus in BaFe2(Asi_2:Pa;)2- The P-for-As substitution results in shrinkage of 
lattice, especially for the FeAs block layers. The resistivity anomaly associated with the 
spin-density-wave (SDW) transition in the undoped compound is gradually suppressed 
by the P doping. Superconductivity with the maximum Tc of 30 K emerges at x=0.32, 
coinciding with a magnetic quantum critical point (QCP) which is evidenced by the 
disappearance of SDW order and the linear temperature-dependent resistivity in the 
normal state. The Tc values were found to decrease with further P doping, and 
no superconductivity was observed down to 2 K for x > 0.77. The appearance of 
superconductivity in the vicinity of QCP hints to the superconductivity mechanism in 
iron-based arsenides. 
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1. Introduction 

The discovery of liigli-temperature superconductivity in iron arsenides [H, E] was an 
exciting breakthrough in the year 2008, [3J which has led to a number of new 
superconductors [H El E] all containing antifluorite-like FeAs or FeSe layers. In most 
cases, superconductivity was induced by doping charge carriers into a stoichiometric 
parent compound which shows a collinear antiferromagnetic spin-density-wave ground 
state[7j. The chemical doping was initially carried out only outside the FeAs 
layers, [Hill [HI [9] with the electronic phase diagram[T0l[TTl[T2] resembling that of cuprate 
superconductors. Later, however, it was found that the even doping at the Fe site (with 
Co [13] or Ni[Tl]) could induce superconductivity. Besides, superconductivity was able 
to be stabihzed in the parent compounds simply by applying hydrostatic pressures. [T5] 
The latter two results contrast sharply with those in cuprates. 

Very recently, we observed superconductivity at 26 K in EuFe2As2 by the partial 
substitution of arsenic with phosphorus. [I6j This isovalent substitution does not 
introduce additional electrons or holes into the system, thus it can be roughly understood 
in terms of chemical pressures. Owing to the ferromagnetic ordering of Eu^^ moments 
below 20 K coexisted with superconductivity, the Meissner effect was not obvious. 
Subsequently, we found that the P/As substitution induced superconductivity in a 
prototype "1111" material LaFeAsO, but the superconducting critical temperature 
Tc is only ~10 K,[T7] which is substantially lower than the T^. value under physical 
pressures [18]. Here we employ the same doping strategy on another prototype "122" 
parent compound BaFe2As2 [19] . Bulk superconductivity up to 30 K with strong 
Meissner diamagnetism has been observed for x=0.32. Our systematic measurements 
of resistivity and magnetization reveal an electronic phase diagram characterized by 
a dome-like superconducting region located near a magnetic quantum critical point 
(QCP). The present work is consistent with a recent theoretical prediction which shows 
a unique type of QCP resulting from a competition between electronic localization and 
itinerancy. [20] 

2. Experimental 

Poly crystalline samples of BaFe2(Asi_a;Px)2 (with the nominal P content a:=0, 0.1, 0.2, 
0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.6, 0.7, 0.8 and 0.9 and 1) were synthesized by solid state 
reaction with the high-purity elements Ba, Fe, As and P. Mixture of the elements in the 
stoichiometric ratio was pressed into pellets in an argon-filled glove-box. The pellets 
were sealed in evacuated quartz tubes and annealed at 600 K for 12 h, then at 1173 
K for 40 h. The resultant was reground, pelletized, and annealed in evacuated quartz 
tubes at 1273 K for 24 h. Finally the samples were cooled down to room temperature 
by shutting off the furnace. 

Powder x-ray diffraction (XRD) was performed at room temperature using 
a D/Max-rA diffractometer with Cu-K^ radiation and a graphite monochromator. 
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The structural refinement was performed using tlie programme RIETAN 2000.pT] 
Tlie electrical resistivity was measured using a standard four-probe metfiod. Tlie 
measurements of dc magnetization were performed down to 1.8 K on a Quantum Design 
Magnetic Property Measurement System (MPMS-5). Both the zero-field-coohng (ZFC) 
and field-coohng (FC) protocols were employed under the magnetic field of 10 Oe. 

3. Results and discussion 

BaFe2(Asi_KPa;)2 system crystallizes in ThCr2Si2-type structure with the space group 
lA/mmm. Except the main phase, minor secondary phase was identified as Fe2P. 
Figure 1(a) shows the XRD pattern and its Rietveld refinement profile of a typical 
BaFe2(Asi_a;Pa:)2 sample with x=0.35 at room temperature. The two-phase Rietveld 
refinement is successful, as indicated by the fairly good reliability factors i?^p=8.5% and 
5*= 1.8. It was found that the fitted P content is a little lower than the nominal one, 
consistent with the existence of impurity Fe2P (about 10% in molar ratio). The refined 
lattice parameters plotted as functions of x are shown in Fig. 1(b). The cell volume 
decreases monotonically with the P doping. We notice that, for the samples with no 
significant Fe2P impurity {e.g., x=0.45 and 0.9), the data points of the cell volume fall 
exactly on the straight line, in accordance with the Vegard's law. Therefore, the actual 
P content can be estimated from the measured cell volume, which gives the corrected P 
content x'=0.02, 0.13, 0.2, 0.24, 0.32, 0.36, 0.47, 0.56, 0.63 and 0.77 for the nominal x 
values of 0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.5, 0.6, 0.7 and 0.8, respectively. This correction 
is basically (within the uncertainty of ±0.02) consistent with the Rietveld refinement 
results. 

Table 1 compares the crystal structures between the undoped and P-doped 
(a;'=0.32) BaFe2As2. A remarkable difference is seen from the position of As/P. The 
P doping makes the pnictogen atoms move toward the iron sheets, which leads to a 
dramatic decrease of FeAs-layer thickness. This explains why the decrease in c-axis 
is more pronounced (compared with the a-axis), as shown in Fig. 1(b). In fact, the 
spacing of FeAs layers decreases little with P doping. This indicates that P doping 
impacts mainly on the FeAs block layers. As a result, the bond angle of As-Fe-As along 
[110] directions increases substantially. 

Figure 2 shows the temperature dependence of resistivity for BaFe2(Asi_i:Pa;)2. For 
the undoped BaFe2As2, p drops rapidly below 140 K, consistent with previous report [19]. 
The resistivity anomaly was revealed to be associated with a simultaneous structural 
and spin-density-wave (SDW) transition. [22] Upon doping with P, the transition is 
suppressed, exhibiting the decrease of transition temperature (Tsdw) and the weakening 
of the resistivity drop. For example, Tsdw is suppressed to ~ 90 K for the sample 
of x'=0.13, where only a resistivity kink appears. When the P content increases to 
20%, a sudden decrease in resistivity is seen below 10 K in addition to an observable 
kink at ~ 70 K. The resistivity drop at lower temperature is due to a superconducting 
transition (confirmed by the magnetic measurements below). In the intermediate doping 
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Figure 1. (a) X-ray powder diffraction pattern and its Rietveld refinement profile 
for the BaFe2(Asi„a;P2:)2 (a;=0.35) sample. The secondary minor phase is identified 
as Fe2P. (b) Relative lattice constants (compared with the undoped compound) as 
functions of nominal P content x. The solid straight line obeys the Vegard's law. 



Table 1. Comparison of room-temperature crystal structures between the undoped 
and P-doped (x'=0.32) BaFe2As2 with the space group lA/mmm. The atomic 
coordinates are as follows: Ba (0, 0, 0); Fe (0.5, 0, 0.25); As/P (0, 0, z). 



Compounds 


BaFe2As2 


BaFe2(Aso.68Po.32)2 


a (A) 


3.9612(3) 


3.9231(3) 


c(A) 


13.006(1) 


12.805(1) 


V (A^) 


204.08(2) 


197.08(2) 


z of As/P 


0.3546(2) 


0.3530(2) 


FeAs-layer thickness (A) 


2.721(2) 


2.638(2) 


FeAs-layer spacing (A) 


3.782(2) 


3.765(2) 


As-Fe-As angle (°) 


111.0(1) 


112.2(1) 
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Figure 2. Temperature dependence of resistivity for BaFe2(Asi_j:'Pa;')2 samples. The 
data are normalized for the convenience of comparison. All the curves are labeled with 
the corrected P content. The insets display expanded plot for showing the fit of the 
normal-state resistivity with the formula p{T) = pQ + aT". 



range (0.24< x <0.63), superconductivity with higher T^s was observed. The maximum 
transition temperature is 30 K for a;'=0.32. The resistive transition tail is probably 
due to the influence of grain boundaries and/or minor ferromagnetic Fe2P impurity, 
rather than sample inhomogeneity, since nearly perfect diamagnetism was measured 
(see below). For the higher P doping (x >0.77), no superconductivity was observed 
down to 3 K. One notes that the residual resistance ratio (RRR), defined by the ratio of 
the resistance at room temperature to the residual resistance at very low temperature, 
is 20, 35 and 75 for x'=0.77, 0.9 and 1.0, respectively. The high RRR values suggest 
high quality of the samples. 

To confirm superconductivity, the temperature dependence of dc magnetic 
susceptibility for BaFe2(Asi_^Pi^)2 was measured, as shown in Fig. 3. The sample 
of x'=0.2 shows strong diamagnetic transition below 7 K, consistent with the drop in 
resistivity below 10 K shown in Fig. 2. In spite of the large superconducting fraction, 
however, it is difficult to ascertain whether SDW order coexists microscopically with 
superconductivity or not. For x=0.32, a stronger diamagnetic signal is observed below 
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Figure 3. Superconducting diamagnetic transitions in BaFe2(Asi-2:'P2:')2. Note that 
the magnetic background coming from the ferromagnetic impurity of Fe2P has been 
removed. No correction of demagnetization has been made. 

30 K. The magnetic shielding is even beyond the ideal value of 100% apparently, since 
the demagnetization effect is not included here. The Meissner volume fraction achieves 
~ 30%, demonstrating the bulk nature of superconductivity. Other superconducting 
samples also displays similar bulk property. 

As is known, the normal-state resistivity contains valuable information on electronic 
scattering and even superconducting mechanisms. Dense polycrystalline samples may 
show intrinsic transport property because the grain boundaries scatter electrons like 
impurities. Thus the following analysis is valid even though our data were obtained from 
polycrystalline samples. Assuming that the intrinsic scattering produces a power-law 
resistivity at low temperatures, the total resistivity can be expressed as p{T) = po + aT"', 
where po represents extrinsic resistivity coming from the impurity/ defect scattering. By 
fitting the p(T) data in the range of Tc < T <50 K, the exponent n was obtained. 
The inset of Fig. 4 shows that most samples have the n-value close to 2.0, suggesting 
Fermi-liquid state due to electron-electron interactions. Interestingly, n is close to unity 
(i.e., linear T-dependent normal-state resistivity) for the samples with a;'=0.3~0.45. 
The linear low-temperature resistivity belongs to a non-Fermi liquid behaviour, as 
usually observed in the regime of quantum criticality. Remember that the SDW order 
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Figure 4. Electronic phase diagram in BaFe2(Asi_a;/P2:')2- SDW and SC denote spin- 
density- wave and superconductivity, respectively. The inset shows the variations of the 
exponent n in the formula p{T) = po + olT^ , obtained by the fitting of normal-state 
resistivity at low temperatures. 

is suppressed with the P doping, we argue that there exists a QCP at x' ~l/3 in 
the BaFe2(Asi_a;Px)2 system. This result is quite consistent with a recent theoretical 
prediction which shows a unique type of QCP resulting from a competition between 
electronic localization and itinerancy[2n]- We aware that similar QCP was also proposed 
in (K.Sr)Fe9As9p!2] system. 

The electronic phase diagram is concluded in Fig. 4. The undoped BaFe2As2 
has the ground state of SDW order. With the P doping, Tsdw decreases gradually, 
forming a phase region of SDW. The suppression of SDW can be easily explained in 
terms of Ji — J2 model[23l El], because J2 is expected to decrease upon P doping. On 
the other hand, the scenario of Fermi surface nesting [25] is not quite straightforward 
to meet the experimental observation. In the vicinity of above-mentioned QCP, an 
asymmetric superconducting (SC) dome emerges. The optimal doping level (with 
the maximum Tc) coincides with the QCP, where a non-Fermi-liquid behviour shows 
up. This strongly suggests that spin fluctuations play an important role for the 
superconductivity. In the overlapped area of SDW and SC, it is not clear whether 
SDW coexists microscopically with SC at this moment. The upper-right metallic region 
displays Fermi-liquid-like state, except that non- Fermi-liquid behaviour appears near the 
QCP. This crossover in electronic excitation states resembles those of (K.Sr)Fe9As9p!2] 
and Ba(Fei_^Co^)9As9[26] systems. 

It was reported that BaFe2As2 became superconducting up to 29 K under the high 
pressure of ~4 GPa. [27] The Tc value is very close to the present one in BaFe2(Asi_a;Pa;)2- 
This implies that chemical pressure generated by P doping acts in a similar way 
to the physical pressure. Meanwhile, the disorder induced by the P doping hardly 
influences superconductivity. The relatively low T^°-^ value (compared with 38 K 
in Bao.6Ko.4Fe2As2[l|) can be well understood by an empirical structural rule for 
Tc variations. [11] A very recent work[28] also reveals similarities between structural 
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distortions under pressure and chemical doping in BaFe2As2. Therefore, it is the 
structural details of FeAs layers (corresponding to specific electronic structures), rather 
than charge doping level, that controls the ground states in iron arsenides. 

To summarize, BaFe2As2 has been tuned superconducting through an isovalent 
doping of phosphorus at the arsenic site, for the first time. Unlike previous doping 
strategy with changing the concentration of charge carrier, the present P doping does 
not alter the number of valence electrons. Structural measurements indicates that 
the chemical pressure generated by the P doping impacts mainly on the FeAs layers. 
The P-for-As substitution suppresses the SDW ordering, and favors superconductivity, 
similarly to phenomena of the carrier doping [1[ [101 HI]- A magnetic QCP is evidenced 
simultaneously by the disappearance of SDW ordering and the non-Fermi-liquid 
behaviour. The maximum Tc of 30 K is observed where the QCP is located. Theoretical 
studies[29] indicate that superconductivity around magnetic QCP is associated with the 
spin fluctuation mechanism. Thus the present result suggests that spin fluctuations play 
crucial role for the superconductivity in iron-based arsenides. 
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